%target
bfp =[4,2, 5, 15,44, 104, 201, 288, 337, 338, 306, 256, 217, 174, 139, 116,97,77,57,45,39,35,27,22,14,11,6,5,3,2,2,2];
mct =[18,6,19,52, 132, 244, 353, 428, 410, 351, 279, 207, 150,92,55,38,26,21,17,13,11,9,10,7,7,6,5,4,5,3,2,5];
bfpPI=[2,  1,  3,  7, 15, 33, 80, 151, 206, 242, 236, 232, 225, 230, 221, 217, 187, 168, 142, 117, 86, 64, 45, 29, 17, 11,  7,  3,  2,  1,  1,  3];
mctPI=[77, 40,85,145,226,302, 346, 332, 304, 255, 225, 179, 134, 100, 71, 47, 39, 30, 14, 13,  8,  4,  2,  3,  1,  0,  1,  0,  0,  0,  0,  0];
bfpEC=[6,  5, 11, 31, 78, 200, 399, 552, 513, 412, 282, 186, 111, 75, 51, 32, 19, 14,  7,  3,  1,  0,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0];
mctEC=[29,  9, 20, 58, 143, 290, 499, 536, 381, 240, 159, 107, 84, 62, 53, 39, 34, 30, 27, 24, 23, 19, 18, 23, 16, 15, 18, 10,  7,  5,  3,  1];

%candidate
%bfpx1 = [12 0 0 0 0 0 0 0 39 0 57 84 99 108 270 264 648 702 525 246 18 0 0 0 0 0 0 0 0 0 0 0 ];
%mctx1 = [228, 0, 0, 0, 0, 0, 0, 0, 294, 0, 282, 261, 264, 279, 438, 336, 477, 138, 66, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
%bfpx2 = [57 0 0 0 0 0 0 0 150 0 231 282 255 249 441 312 522 285 195 87 6 0 0 0 0 0 0 0 0 0 0 0 ];
%mctx2 = [36, 0, 0, 0, 0, 0, 0, 0, 90, 0, 132, 132, 162, 204, 453, 462, 795, 423, 168, 15, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
%bfpx3 = [6 0 0 0 0 0 0 0 21 0 57 81 81 99 252 249 699 594 660 240 33 0 0 0 0 0 0 0 0 0 0 0 ];
%mctx3 = [210, 0, 0, 0, 0, 0, 0, 0, 282, 0, 273, 303, 279, 243, 468, 375, 408, 171, 60, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];

%fpx1 = [15 0 0 0 0 0 0 0 45 0 72 87 102 93 165 147 267 216 309 480 531 396 141 6 0 0 0 0 0 0 0 0 ];
%ctx1 = [384, 0, 0, 0, 0, 0, 0, 0, 432, 0, 420, 363, 267, 234, 378, 246, 201, 129, 15, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
%fpx2 = [99 0 0 0 0 0 0 0 222 0 312 291 249 231 288 225 303 153 165 210 177 114 33 0 0 0 0 0 0 0 0 0 ];
%ctx2 = [63, 0, 0, 0, 0, 0, 0, 0, 69, 0, 120, 87, 72, 96, 213, 219, 462, 570, 708, 324, 66, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
%fpx3 = [9 0 0 0 0 0 0 0 12 0 30 39 24 45 51 54 78 1602 132 159 255 336 228 18 0 0 0 0 0 0 0 0 ];
%mctx3 = [201, 0, 0, 0, 0, 0, 0, 0, 252, 0, 228, 171, 120, 1644, 186, 102, 132, 27, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];

bfpx1 = [0 6 9 54 99 171 243 225 246 237 222 237 237 366 378 282 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 48 ];
mctx1 = [0, 0, 9, 42, 123, 198, 240, 249, 252, 204, 228, 222, 294, 333, 363, 261, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 48, ];

bfpx1 = [0 0 6 6 51 132 255 300 288 297 369 330 435 366 192 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 36 ];
mctx1 = [0, 0, 3, 15, 60, 141, 264, 264, 339, 345, 330, 339, 345, 423, 150, 18, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 36, ];

bfpx1 = [0 3 0 15 51 159 249 327 309 282 366 351 411 357 144 18 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 30 ];
mctx1 = [0, 3, 3, 12, 57, 138, 228, 300, 297, 330, 309, 399, 381, 405, 177, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 30, ];

bfpx1 = [0 0 0 3 90 132 219 252 273 294 339 321 474 420 207 15 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 33 ];
mctx1 = [0, 0, 3, 27, 66, 192, 237, 345, 318, 312, 351, 306, 354, 339, 180, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 33, ];

bfpx1 = [0 0 0 12 57 144 291 375 285 225 198 264 444 612 126 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 39 ];
mctx1 = [0, 0, 3, 6, 48, 141, 318, 426, 285, 234, 177, 270, 447, 531, 141, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 39, ];

bfpx1 = [0 0 0 6 45 165 267 339 339 246 324 417 612 246 18 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 48 ];
mctx1 = [0, 0, 3, 6, 30, 138, 261, 471, 354, 249, 291, 366, 543, 291, 21, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 48, ];

bfpx1 = [0 0 3 18 75 198 273 321 363 369 261 255 228 276 177 117 81 12 3 0 0 0 0 0 0 0 0 0 0 0 0 42 ];
mctx1 = [0, 3, 0, 33, 117, 192, 306, 372, 360, 360, 276, 237, 177, 204, 231, 111, 42, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 42, ];
figure(1)
clf

subplot(2,1,1); hold on;
title('BFP (Uninduced)')
plot(bfp,'b')
plot(bfpx1,'r')

subplot(2,1,2); hold on;
title('mCitrine (Uninduced)')
plot(mct,'b')
plot(mctx1,'r')
